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We study fermionic superfluidity in a boson-single-species-fermion cold atom mixture. We argue 
that apart from the standard p-wave fermion pairing mediated by the phonon field of the boson 
gas, the system also exhibits s-wave pairing with the anomalous correlator being an odd function 
of time or frequency. We show that such a superfluid phase can have a much higher transition 
temperature than the p-wave and may exist for sufficiently strong couplings between fermions and 
bosons. These conditions for odd-frequency pairing are favorable close to the value of the coupling 
at which the mixture phase-separates. We evaluate the critical temperatures for this system and 
discuss the experimental realization of this superfluid in ultracold atomic gases. 
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I. INTRODUCTION 

Symmetry of superfluid phases in correlated fermion 
liquids has been a subject of extensive research for many 
decades. Two of the most famous examples of systems 
exhibiting non-trivial symmetries of the order parameter 
are superfluid 3 Hc and high-T c cuprates, where pairing 
is believed to occur in p-wave triplet and d-wave sin- 
glet channels, respectively. In 1974, in an attempt to 
explain the existence of several superfluid phases in 3 He, 
Berezinskiii had suggested that there is yet another pos- 
sibility for the triplet pairing: while the s-wave com- 
ponent of the anomalous correlator (ip a ( r : t)^p{ v> 1 1')) 
should identically vanish at equal times due to the Pauli 
principle^., it may still be nonzero for t ^ t' (i.e., the 
anomalous correlator is an odd function of time or fre- 
quency) thus giving rise to the superfluid 3 He A-phase. 
While experiments have shown that Berezinskii's conjec- 
ture is not realized in the actual 3 He liquid, the idea of 
such nonlocal-in-time pairing has attracted considerable 
attention later, with the discovery of superconductivity in 
cuprates 3 ^ and heavy fermion compounds^. In particu- 
lar it has been argued that the nonlocal character of such 
an order parameter provides a natural resolution to the 
"paradox" of the coexistence between the strong, but in- 
stantaneous short-range Coulomb repulsion and Cooper 
pairing in these materials. More recently the odd pairing 
mechanism was proposed to explain the anomalous prox- 
imity effect in superconductor-fcrromagnet junctions?, 
where superconducting penetration length is believed to 
be significantly enhanced due to the formation of the 
triplet odd-parity component from the standard s-wave 
singlet condensate. 

Despite considerable interest in the subject, the ques- 
tion of whether an odd-frequency phase exists in equilib- 
rium physical systems with no external source of pairs, 
i.e., as a consequence of spontaneous symmetry break- 
ing, remains unresolved. Previous work has suggested 
that such phases may be thermodynamically unstable^, 
but the situation remains unclear. Here we discuss the 
possibility that a boson-fermion mixture, presently real- 
izable in atomic traps, provides an example of a system 



where an odd-frequency fermionic superfluid phase may 
exist under the appropriate conditions. 

We show that due to the interaction with the phonons 
in the bosonic subsystem, the fermions at sufficiently low 
temperatures exhibit pairing either in the p-wave channel 
or in the s-wave odd-frequency channel. A key result is 
that the s-wave odd-frequency pairing exists only when 
the coupling 7 (to be defined below) between the phonons 
and the fermions exceeds a certain threshold value j c . 
Moreover the value of 7 C is close to the coupling strength 
at which the mixture phase-separates. That is, upon 
an increase in the boson-fermion coupling the phonon 
mode softens, thus leading to stronger attractive interac- 
tion between the fermions, as a result of which the odd- 
frequency fermionic condensate can form in the vicinity 
of the phase separation transition. We estimate the tran- 
sition temperatures for the system. We also point out 
certain cross-correlations between the boson and fermion 
densities that could possibly be detected in cold atom 
mixtures as a signature of the odd-frequency superfluid 
phase. 



II. MODEL 

We describe the dilute mixture of fermions and bosons 
with the following Hamiltonian density: 



H = H° B + ^y B \ 



(1) 



where H B F denote the Hamiltonians for noninteracting 
bosons and fermions, A and A' are the boson-boson and 
boson-fermion coupling constants. We assume that the 
trap confining the particles is magnetic and therefore the 
fermions are fully spin polarized. As a result, due to the 
Pauli principle, direct interaction between fermion atoms 
is negligible. For the purposes of the present calcula- 
tion we also neglect the spatial variation of the trapping 
potential and assume that the local fermion and boson 
densities are controlled by the chemical potentials so that 
H BF = i/j b F (— V 2 /2m_B,_F — [<tB ,f)^1'b ,f ■ Also, here and 
in the following we set h = k B = 1. 
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In order to derive an effective coupling between 
fcrmions it is convenient to rewrite the bosonic fields in 
terms of the density-phase variables as Vb = y^exp (i(j>). 
Furthermore, writing p = po + dp, where po is a con- 
stant and Sp contains only nonzero frequency tu n and 
wavevector q components of p{uj n , q), expanding Eq. ([1]) 
up to 0(5p 2 , 4> 2 ) around po, and integrating out the phase 
variable, we obtain the effective "electron" -phonon model 
with the Matsubara action 

S cS = J" ^ dr J d 3 r (C% + X'Sp^Uf) (2) 

2 V p <r 4m B poJ 

In Eq. © C F is free-fermion Lagrangian, j3 = 1/T 
where T is the temperature, and the last term describes 
phonons or Bogolubov quasiparticles with dispersion re- 
lation ui = c s q(l + Co*? 2 ) 1 ^ 2 ' where the phonon speed 
of sound c s = (Xpo/mB) 1 ' 2 and the boson coherence 
length £o = (l/4mBApo) 1//2 - The phonons give rise to the 
nonlocal fcrmion density-density interaction, with vertex 
— X D°(u! n , q)/2, where D 0_i (w„,q) is the expression m 
the brackets in the last term in Eq. @. 

It should be noted that while the effective description 
in Eq. @ is valid for finite values of the coupling con- 
stant A', it breaks down for sufficiently large A' or fermion 
density. Indeed, by evaluating the rcnormalization of the 
phonon Green's function D° within the second order (in 
A') perturbation theory in the particle- hole channel (i.e., 
accounting for a single fermion polarization bubble in 
phonon self-energy), we see that 

D - 1 D" 1 = ^ + A(l - 7) + A£V + 0(q% (3) 

Pot 

where 7 = X' 2 q 2 F /(2TT 2 Xv F ) and £ 2 = $ + 7/12q>|; q F 
and v F = qF/m-F are the Fermi momentum and velocity 
respectively^. Thus for 7 — > 1 the phonon mode softens 
at small q, signaling that the mixture becomes unsta- 
ble against the phase separation transition. The value of 
7 = 1 corresponds to the line of spinodal decomposition 
in the mean field analysis where the instability shows up 
as a saddle point in the free energy-^. While the second 
order perturbation theory will not be quantitatively accu- 
rate for 7 close to 1, the higher order renormalizations of 
the phonon propagator do not alter this conclusion on the 
qualitative level; their effect merely leads to a redefini- 
tion of 7 and As a result effective interaction between 
fermions increases as one approaches the spinodal point, 
diverging for u> = q = at 7 = 1. It should be pointed 
out that the phase separation is, presumably, a first order 
phase transitio n 10 ' 12 , and therefore the separated phase 
becomes thermodynamically more favorable before the 
spinodal line. However, an estimate of nuclcation rates 
shows that due to extremely low temperatures as well as 
relatively weak intcrparticle interactions such nuclcation 



processes are exponentially slow. That is, upon an adi- 
abatic increase of fermion-boson interaction the system 
will remain in the metastable mixed state up until the 
absolute instability (spinodal) line, unless the boson gas 
parameter is comparable to 1, see RefJ^. 



III. GAP EQUATION 

The onset of the pairing instability corresponds to the 
appearance of a nonzero anomalous correlator F(r — 
r',r — t') = i(Tip F (r,T)ip F (r' ,r')). The self-consistency 
equation can be readily obtained within the Eliashbcrg 
formalism^. Following the standard procedure we derive 
the linearized eigenvalue-type equation for F(u>, q): 

G-'K.qjG-'U.qjFK.q)^ £ F(oj' n ,q') (4) 

A' 2 

x — [D("n - u'„, q - q') - D(v n + u' n , q + q')] > 

where G is the fermion Green's function and u> n = 
nT(2n + 1) is the Matsubara frequency. Note that the 
appearance of the difference in the brackets on the right 
hand side (r.h.s.) of Eq. (fj| is due to the presence of 
fermions with the same spin only; had we considered the 
usual singlet pairing, the difference would have been re- 
moved by the spin part of the Cooper pair wavefunction. 
As a result solutions to Eq. ^ should satisfy the anti- 
symmetry property, F(uj n , q) = — F{— ui n , — q), which, as 
expected, rules out a possibility of the standard s-wave 
even-in-time pairing. 

In order to solve Eq. ^ we must specify the Green's 
functions G and D, which are renormalized in the 
particle-hole channel according to the Dyson equations 
G- 1 = G 0_1 - E and D^ 1 = D ^ 1 - II, where E and IT 
are the fermion and phonon self-energies, e.g., Eq. ([3]). 
It is well known from the Fermi liquid theory^ that 
near the Fermi surface the renormalized fermion Green's 
function G can be written as Z[iuj n — v F (q — qp)]^ 1 , 
where Z is the quasiparticle "weight" coefficient [Z < 
1), and v* F = Zvp. Then G _1 (w„, q)G _1 (— u> n , q) = 
[uj 2 + v* F 2 (q - q F ) 2 ]/Z 2 . Next let us expand F(w,q) 
and D(cu,q) in orbital harmonics using Legendre poly- 
nomials Pi . For reasons to be specified below we assume 
that c s /v F is small. Then we notice that D is a rela- 
tively slowly varying function of q compared to F (which 
is strongly peaked at q = q F ) and therefore we can set 
both |q| and |q'| in the D's on the r.h.s. of Eq. (|2|) equal 
to q F . As a result the Cs are functions of the angle be- 
tween q and q' only, and with the use of the addition 
theorem we obtain 

K + v f ]Fl(wn,8q) = — - — T 2^ J (5) 
x^K, 5q') [T>i(w n -u' n ) - (-l)%(LJ n + u>' n )] , 
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where Ti is the partial component of the anomalous cor- 
relator F (i.e., with angular momentum I), Sq = q — qp, 

and Vi = j\ dcos6D(oj,q F y/2 - 2 cos 6) Pi (cos 6). 

It is obvious that there are two types of solutions to 
Eq. ([5]): !Fi{lo) with odd I and even in to and vice versa, 
with even I and odd in ui. For both of these solutions the 
bracket in the r.h.s. of Eq. ((SJ) can be replaced by 22?; (w— 
u)') and then we note that Eq. ((SJ can be cast in the 
form HiTi = 0, where Hi is a "Hamiltonian" of a particle 
moving in a two-dimensional external potential Vi (r, a;) ~ 
— 2T>i(t)6(x) (here r is Matsubara time and x is Fourier 
conjugate to Sq), and J r ;(/3/2, x) = JF/(— /3/2, x) = 0. 

For /3 — ► oo, i?/ has at least one negative eigenvalue 
£oo,z corresponding to the bound state of the particle. 
Upon an increase of T (i.e., for finite 0) the "energy" 
£oo,i increases, reaching for some particular value of 
(3 = (3 C . In order to estimate the value of j3 c we note 
that for compact "Di(t) the bound state "wavefunction" 
behaves asymptotically as ~ exp(— e^r), and there- 
fore the boundary conditions, i.e., the finitencss of /3, 
start to intervene when (3 ~ 1/-^/— £oo,z. Moreover, since 
(— £c»,;) _1 ^ 2 i s the only relevant "time" scale (in the limit 
when (3~ l is much smaller then the cut-off frequency oj c 
of T>i; u> c ~ c s qp), the critical temperature T Ci z = l//3 c ,z 
for channel ^ is of the order of -J—t^i- 

Evidently the ground states of fl/'s describe only the 
channels with odd I and even w-dcpcndcncc of J-iS. The 
second type of solutions, even in I and odd in r, how- 
ever, correspond to the first excited states of Hi, which, 
obviously, are odd functions of uj or r. While the same 
considerations hold for the critical temperatures for these 
solutions as well, there is an important distinction: in the 
(3 — > oo limit the former solutions (e.g., bound ground 
states) always exist for the attractive Vi's, the latter 
may only exist when the potential Vi is strong enough. 
Therefore even at zero temperature the odd-frequency 
(Berezinskii) states exist only if the coupling constant A' 
or 7 exceeds a certain threshold j c . Moreover, since at 
the threshold point the state's "wavefunction" is delocal- 
ized and it localizes for stronger Vi, the odd- frequency 
phase emerges at zero temperature when 7 = 7 C and ex- 
tends to non-zero temperatures at stronger couplings. 

It is easy to estimate this threshold coupling 7 C and 
the critical temperatures from Eqn. ([3]); see appendix. 
For the I = odd frequency case, the zero temperature 
critical coupling 7 C can be estimated as the solution to 
7c = (a 2 /Z 2 )/\n[l + a 2 /(l - 7c )], where a = 2q F £; 
see appendix. The solutions to this equation arc ly- 
ing in the interval 1/(1 + Z 2 ) < j c < 1, depending 
on the value of a. In particular, for a 3> 1 one finds 
7 C s=! 1 — a 2 cxp(—a 2 /Z 2 ), and so the width of the super- 
fluid phase (in terms of the dimcnsionless coupling 7) is 
exponentially small. For a ~ 1, 1 — 7 C is finite and thus 
the odd-superfluid phase is most favorable when the ef- 
fective boson coherence length £ is comparable with the 
interfermion distance. The transition temperature can 
be calculated numerically, sec appendix, and the results 
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FIG. 1: (a) thru (c): Critical temperatures vs. fermion- 
phonon coupling 7 for a — 2qp^ = 1.0, 1.4, and 2.0. The 
solid line is the critical temperature for s-wave odd-frequency 
pairing and the dashed line for p-wave pairing. In all cases we 
have set Z = 1. The z-axis is the same for all three panels, 
but note the different j/-axis scales. 



are shown in Fig. l(a-c) for different values of a. Note 
that for a ~ 1 the critical temperature reaches values 
~ c s qF- 

For the even frequency p-wave case, we consider the 
effective potential in the I = 1 channel, using T>\. We 
note that this is the same as the p-wave case studied in 
RefsJ^ except that we use the renormalized T>\ following 
from Eq. ([3]) — this places the p-wave on equal footing 
for comparison to the the s-wave odd frequency phase. 
The p-wave phase exists as the coupling 7 goes to zero, 
although the critical temperature is exponentially sup- 
pressed. The p-wave critical temperature is plotted in 
Fig. l(a-c) for different values of the parameter a as a 
dashed line. Note that due to the relative smallness of 
the effective coupling strength T>\ (0) the critical temper- 
ature for the p-wave pairing is much lower than that for 
the odd frequency s-wave (except for Fig. 1(c), where the 
odd- frequency region is very small, i.e., for large a). 



IV. POSSIBLE EXPERIMENTAL 
REALIZATION AND DISCUSSION 



Since most of the standard (in solid state) thermo- 
dynamic and transport measurements are presently un- 
available in cold atoms, detecting odd-frequency pairing 
in a cold atom fermion-boson mixture may be a non- 
trivial problem. One intriguing possibility would be to 
take advantage of the cold atom time-of-flight type ex- 
periments to study the unique correlations of the odd- 
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FIG. 2: Schematic picture of three-particle correlations under 
time-of-flight expansion: two fermions, rw,i and rF,2, and one 
boson, tb. 



frequency phase. For example, it has been demonstrated 
in Ref that momentum correlations of atomic fermions 
can be observed by the photodissociation of molecules; 
upon release from a trap, the atoms exhibit density cor- 
relations between points r and — r (relative to the center 
of the trap) as a consequence of the initial molecular 
state. Similar measurements have been proposed to de- 
tect other types of many-body correlations^; for exam- 
ple, a fermionic gas in the BCS regime has a density- 
density correlation function (np(ri, t)riF(^2, t)) which 
at sufficiently large times of flight t is also peaked at 
ri = — Y2 as a result of the pairing. 

The above argument should, however, be modified 
when applied to the odd-frequency type of pairing. In- 
deed, since the equal time anomalous correlator is iden- 
tically for the odd-frequency pairing, there is no 
V\ = —V2 correlation in the two-point fermion correlation 
function. However, the three-point correlation function 
(ni?(ri)7ii?(r2)ns(r3)), where is boson density, does 
contain the signature of the odd-frequency pairing. To 
see this one should notice that while (V'Fq(O)V'F-q(O)} = 
for such pairing, (^Fq(0)V>F-q(0)) ^ 0. Then, using 
iipp — [H, V>f], where H is given by Eq. (JD), it is easy to 
show that (V'Fq 1 (O)-0Fq 2 (O)ri Bq3 (O)) ~ <5(qi + q 2 + q 3 ), 
where nsq is Fourier component of boson density and 
the width of 6 is primarily controlled by the size of the 
trap. The boson density can be expressed as fifl q « 

y / po(V'Sq +V'B-q)j an< ^ therefore the three point correla- 
tion function contains an irreducible contribution peaked 
at qi + q2 = — q3 (in Fourier space). This three-point 
correlation can be interpreted as an order parameter first 
proposed in the context of odd-frequency superconduc- 
tivity in a t-J modeU&. As a result the real space equal- 
time correlation function (np{vi)np{v2)nB has a 
correlation peak at ri + r 2 = — (m^/fiF)^, where we 
have used the relationship between the wavevectors of 
the particles in the initial state with their coordinates in 
the time-of-flight image, qF,s = ^F.sr/i; see Fig. 2. 
Therefore particle density cross-correlations which can, 
in principle, be deduced from instantaneous fermion and 
boson atom absorption images would provide a direct test 
for odd-frequency pairing in cold-atom mixtures. 



Finally we point out that the validity of Eq. ([4]) is 
controlled by the Migdal criterion, i.e., the smallness of 
the vertex correction part. The latter can be estimated 
as the three-legged diagram 8 -, with two fermion and one 
phonon lines. The standard order of magnitude estimate 
for the ratio between the bare and the renormalized in- 
teraction vertices gives (7/Vl — 7) x (As/Vf) x (^fC)" 1 , 
where we have used the phonon Green's function given by 
Eq. Q. Thus, Eq. Q is valid for finite 7, provided that 
7 is not too close to 1, i.e., the point of the phase sepa- 
ration, and (c s /v F ) x (<7f£) _1 ^ !• Moreover, since we 
are interested in the regime <?f£ ~ 1, the vertex correc- 
tion is small as long as c s /vp <C 1. Therefore the above 
results are quantitatively valid when £oc s rriF <C 1 (K), 
that is mF/mB <C 1. While this condition is naturally 
fulfilled in solid-state systems, with mp and ms being 
electron and ion masses, this is not necessarily the case 
in cold atom systems. A reasonable choice for testing 
our theory in trapped cold atom systems is a 6 Li- 87 Rb bi- 
nary mixture so that the mass ratio condition is satisfied. 
Using the 87 Rb background scattering length of 5.32nm 
and a 6 Li density of 10 12 cm -3 , the coupling 7 = 1 cor- 
responds to a 87 Rb- 6 Li scattering length of 22.7nm, an 
order of magnitude accessible via an interspecies Fesh- 
bach resonance. To optimize the odd-frequency phase, 
one would want a = 2qp£ ~ 1 (see Fig. 1). With the 
same parameters as above, 2qp^a = I corresponds to a 
87 Rb boson density of 2.3 x 10 14 cm -3 and the tempera- 
ture scale T c ~ c s qp ~ 85nK, readily accessible in cold 
atomic gases. 

In summary, we have studied a mixture of bosons and 
singlc-spccics-fcrmions, showing that fermionic superflu- 
idity of the Bcrczinskii odd-frequency type is likely to ex- 
ist under appropriate conditions, i.e., this pairing occurs 
with a finite critical strength of the boson-fermion cou- 
pling. We have estimated the transition temperature and 
pointed out the unique boson-fermion cross-correlations 
which such a state exhibits. 
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VI. APPENDIX 

To evaluate the threshold coupling and the crit- 
ical temperatures it is convenient to transform the 
"Schrodingcr" equation H1T1 = t x ^Ti back to the fre- 
quency representation, e.g., Eq. (5), but with the summa- 
tion over the Matsubara frequencies replaced by integra- 
tion and with an additional £00,1^1 term on the l.h.s. In- 
troducing A t (u), Sq) = [w* + v* F 2 Sq 2 - e 00 ,;]^ 7 i(o; n , Sq), in- 
tegrating out the momentum 8q' on the r.h.s., and notic- 
ing that the solution A; is independent of Sq (a result of 
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evaluating the phonon propagator on the Fermi surface) , 
we obtain 

We consider the I — and I — 1 solutions to Eqn. (IA1I) . 
Since the effective coupling strength P;(0) rapidly de- 
creases for greater I, phases with higher orbital momen- 
tum of the order parameter have much lower critical tem- 
peratures and thus are never realized. 

To estimate the critical coupling strength for the / = 
odd-frequency solution let us first consider the zero tem- 
perature case (eoofi = in Eq. (|A1[0 . Since A (w) 
must be odd, we use the ansatz Aq(uj) oc u> (with cut- 
off cu c ~ c s qp); we have verified that this linear in- 
dependence as u> — » is correct by means of explicit 
numerical solutions of the gap equation. Expanding the 
r.h.s. of Eq. (|A1|) linearly in w, canceling uj and inte- 
grating the resulting r.h.s. of Eq. (|A1|) by parts we ob- 
tain the condition for the threshold (critical) coupling: 
4ir 2 Vp ~ (X' c Zq F ) 2 T> (0) (here we have neglected the 
cut-off dependence assuming that Vo{uj c ) *C X>o(0)). It 
is instructive to evaluate the critical coupling strength 
for -D(O) = D°(0), that is, without accounting for the 
phonon-mode softening. A straightforward calculation 
yields 7 C = \' c 2 q 2 F / {2tt 2 \v f ) = a§/[Z 2 ln(l + a 2 ,)], where 
Q'o = 2gi?£o- Since this "f c is always greater than 1 
(note that Z < 1), one would conclude that the cou- 
pling needed for the formation of the Berezinskii phase 
is stronger than that of the phase separation (7=1) 
and, therefore, that the phase does not exist. This con- 
clusion, however, is erroneous because at finite fermion- 
phonon coupling the renormalization of the phonon prop- 
agator, e.g., Eq. (3), is crucial: as the coupling strength 



approaches that of the phase separation, the effective in- 
teraction between fermions increases due to phonon soft- 
ening. Thus, with D(0) given by Eq. (3) we find that 7 C 
satisfies the equation a 2 /(Z 2 j c ) = In [1 + a 2 /(l — 7 C )] 
where a = 2q F ^, which has solutions 7 C < 1. In particu- 
lar, for a 3> 1 one finds j c ~ 1 — a 2 cxp(— a 2 /Z 2 ), and so 
the width of the supcrfluid phase (in terms of the dimcn- 
sionlcss coupling 7) is exponentially small. For a ~ 1, 
1 — 7 C is finite and thus the supcrfluid phase exists only 
when the effective boson coherence length £ is compara- 
ble with the interfermion distance. The transition tem- 
perature can be estimated within the same linear ansatz 
for Ao(w) by retaining the — £oo,o hi the denominator on 
the r.h.s. and solving the resulting equation numerically 
for yj— £00,0 ~ T c fi, for which the results are presented in 
Fig. l(a-c) for different values of a. 

For the I = 1 even-frequency p-wave phase, to esti- 
mate the critical temperature T C; i or e Mi i from Eq. (|Aip 
it is sufficient to set Ai(w) = const (again, with cut- 
off ~ c s qp) and replace T>\{uj — u>') in the r.h.s. of 
Eq. (|Alj) by 2?i(0). After a straightforward calculation 
using the phonon propagator of Eq. (3) we obtain that 
T c ,\ ~ c s q F exp 

where, for simplicity, we have set Z = 1. We note that 
the p-wave phase has been considered in Refs. [14] for 
fcrmion-boson mixtures, and that our calculation would 
have been the same as those if we had not used the renor- 
malized phonon propagator. The p-wave transition tem- 
perature given by Eqn. (|A2[) is plotted in Fig. l(a-c) for 
different values of the parameter a as a dashed line. 
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